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Abstract. We investigate the order parameter of Sr2Ru04 in an exchange-only mag- 
netic field. A Ginzburg-Landau symmetry analysis implies three possibilities: a pure 
^He A phase, a "^He Ai or a '^He A2 phase. We explore the exchange field dependence of 
the order parameter and energy gap in a one-band model of Sr2Ru04. The numerical 
solutions show no Ai phase and that the A2 phase is lower in free energy than the A 
phase. We explore heat capacity as a function of temperature and field strength and 
find quantitatively different behaviors for the A and A2 phases. 



1 Introduction 

A spin triplet superconductor should show a number of interesting magnetic-field 
effects which are direct consequences of the magnetic moment of the Cooper 
pairs. In particular, for spin-triplet superconductors the Zeeman coupling be- 
tween the quasiparticle spins an external magnetic field need not lead to Pauli 
limiting, unlike the case of spin-singlet superconductors. In the extreme high field 
limit with completely exchange-split bands we could expect single spin pairing 
of the majority spin Fermi surface. We may also expect possible phase transi- 
tions or symmetry changes of the order parameter in a magnetic field, which 
are analogous to the transitions seen in superfluid ^He^H]. The ^He B-phase 
is destroyed in a magnetic field in a qualitatively similar manner to a singlet 
superconductor. On the other hand, if the zero field ground state is one of equal- 
spin-pairing (ESP), then the gap function can deform continuously as a function 
of a Zeeman field the '^He A-phase evolves first smoothly into A2 phase and 
then, via a phase transition Ai phase, as it progresses from equal-spin-pairing 
to single spin pairing with increasing field Q. 

The superconductor Sr2Ru04 should be an ideal candidate to examine these 
effects. There is strong evidence for spin-triplet pairingp) from direct measure- 
ments of spin-susceptibility in the superconducting state|^,^. It is has a simple 
and well understood Fermi surface and is in the clean limit. The detailed 
gap function is still somewhat controversial, but is generally believedj|] to be 
of tetragonal symmetry [p^, and more specifically to be a two-dimensional 
analogue of the '^He A-phase, with d(k) (sin + i sin fcj,)(0, 0, 1). This order 
parameter would agree with the spin-susceptibility measurements in the super- 
conducting state and also would lead to time- reversal symmetry breaking 
below Tc [pdl. More recently, specific heat[p2|, penetration depth |13] and ther- 



mal conductivity! 14 experiments have shown that the gap must have line- nodes 
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on the Fermi surface. However, for the cylindrical Fermi surface geometry |Q 
of Sr2Ru04 a complete group theoretic analysis of symmetry distinct pairing 
states does not show any which both break time-reversal symmetry and have 
line-nodes A possible resolution to this dilemma has been developed in the 
orbital dependent pairing model of Zhitomirsky and Ricejl^ and a in related 
model by Litak et aZ.|l^]. For different reasons both groups proposed that the 
gap function is of the form 

d(k) - (sinfc:r + isinfcj,)(0,0, 1) (1) 

on the dominant 7-Fermi surface sheet, and of the form 

d(k) ^ (sin cos +zsin cos cos (^) (0,0,1) (2) 

on the a and /3 sheets. Both of these functions possess the same symmetry, 
but correspond to intra-plane and inter-plane pairing interactions respectively. 
This gap function has horizontal line nodes at = ±7r/c on the a and (3 
sheets, and was shown to be in good agreement with experimental temperature 
dependences for specific heat, penetration depth and thermal conductivity [p^ . 

In a magnetic field Sr2Ru04 shows a number of unusual features. Firstly 
the vortex lattice is square which agrees well with the predictions of a 

two-component symmetry Ginzburg-Landau theory Secondly there is 
an anomalous second feature close to Hc2, which only occurs when the field is 
aligned within 1° of the a-b plane p^. At the present time the origin of this 
feature is uncertain. It may be a vortex lattice phase transition, or it may cor- 
respond to a change in pairing symmetry with field, perhaps analogous to the 
double superconducting transition in UPts . 

In this paper we will focus specifically on the unique effects of the Cooper 
pair spin in a triplet superconductor. Therefore we neglect the effects of the 
vector potential on the quasiparticles, and instead focus solely on the Zeeman 
coupling of the quasiparticle spin to the magnetic field. We can justify this model 
by appeahng to the strong Stoner enhancement in Sr2Ru04 ||2^ , ^ , and so the 
exchange field will be large. Alternatively, our model may be appropriate for the 
ferromagnetic superconductor ZrZn2 [p4[ . 

This paper is organised as follows. Firstly we write a simple single-band 
model Hamiltonian for p-wave pairing in the 7-band of Sr2Ru04 . Next we 
examine how a spin-only magnetic field enters the corresponding E^ symmetry 
Ginzburg-Landau theory. In section 4, we present detailed numerical results for 
the field dependent energy gap, and specific heat for the two relevant cases of the 
exchange field either parallel or perpendicular to the d(k) order parameter. We 
show that the lower free energy state in analogous to the '^He A2 phase. Finally, 
in Sec. 5 we present our conclusions. 
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2 A Microscopic Model for a Triplet Superconductor in a 
Spin Only Magnetic Field 

We consider the effect of a spin-only magnetic field, H, on an attractive, nearest 
neighbour, Hubbard model. Wc use a one-band model, appropriate for the 7 
sheet of the Sr2Ru04 Fermi surface. The set of interaction constants, U[j'' , 
describe attractions between electrons on sites i and j with spins a and a'. The 
Hamiltonian for this model is: 



where tij is the hopping integral, e is the site energy, c^J are the usual an- 
nihilation (creation) operators and rijcr is the number operator. tTcrcr' are the 
components of the vector of Pauli matrices: 

g. = ^) (4) 

By making the Hartree-Fock-Gorkov approximation and taking a lattice Fourier 
transform the following spin-generalised Bogoliubov-de Gennes (BdG) equation 
can be derived from the above Hamiltonian. 



ijaa' iaa' 

(3) 
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,(5) 



where efc is the (Fourier transformed) normal, spin independent part of the 
Hamiltonian. The order parameters Z\ct<t' (k) are determined self consistently by 



A„,,{\,) = - - ^t/--'(q)(«,,„(-q)<,^„(-q) - v;^„{q)u,,,„{q)){l - 2^.) 

qrr" 

(6) 

where U'^'^' (q) is the lattice Fourier transform of U[j'^' and fqa is shorthand for 
the Fermi function f{Eq„). It is natural to separate the spin-generalized BdG 
equation into triplet and singlet parts: 
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do{k) is the (scalar) singlet order parameter and d{k) is the (vector) triplet 
order parameter. The singlet order parameter is symmetric under spatial inver- 
sion while the triplet order parameter in anti-symmetric under spatial inversion. 
Hence, the BdG equation can be rewritten as 



/ efc + 
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(8) 



If there is no superconductivity in the triplet channel we regain the standard 
result [|j for the spectrum of a singlet superconductor in a spin only magnetic 
field: 



Eik)^±^el + \doik)\^±^^B\H\. (9) 

By setting the singlet order parameter to zero we find that the equivalent result 
for a triplet superconductor is 



Eik) = ± V4 + t^l\H\' + |d(fc)P ± y/m (10) 

where 

A{k) = \d{k) X <i(fe)*|2 + ieltil\H\^ + 4/4|if • d{k)\^ 

+4iekfJ-BH ■ d{k) x d{k)* . (11) 

It should be noted that this does not assume a unitary order parameter]^. 

It is a relatively straightforward process to calculate thermodynamic prop- 
erties for a triplet superconductor. For example the specific heat is given by 

^ A unitary state is any state for which d{k) x d*(k) — 0. 
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Cv = T— (12) 

- -ksT^ Y^ifkaHfka) + (1 - fka)ln{l - fka)) (13) 

= E Ml - A.) (Iff + j^E.(^)^i=»(*)) (") 



3 Ginzburg— Landau Theory of a Quasi— Two Dimensional 
Triplet Superconductor in a Magnetic Field 

Before considering numerical solutions of the self consistent Bogoliubov-de Gennes 
equations, we will examine the possible results by deriving a Ginzburg-Landau 
theory from our microscopic theory. 

Consider a quasi-two dimensional system with two orbital degrees of freedom 
(which we label x and y) and three spin degrees of freedom (labelled 1, 2 and 
3.) Hence, instead of the familiar 3 by 3 order parameter of '^He this system is 
described by the complex 2 by 3 matrix A, which is related to the microscopic 
order parameter, d(fc), by 



sin 



^di{k)\ / Aix sin kx + Aiy sin ky 

d2{k) = A2x sin kx + A2y sin ky j = (A^ Ay) f ) • (16) 

yd3{k)J \A3x sin kx + A3y sin ky 

In zero field the free energy for a tetragonal crystal is given by||lO| 



F^a{T- T,){\Axf + \Ayf) + PiQA^f + \Ay\^f + /^al^. ■ A^ + Ay ■ Ay\^ 
+(3s{{Ax ■ A;r + (A: ■ Ayf + (A* . A,f + (A; • Ayf) 
+/34(2|A,.A;|2 + |A,|4 + |A,|4) 
+/35(2|A, • Ay\'' + \Ax ■ + \Ay ■ Ay\'') 

+ /56(|A, -^,12 + \Ay ■ + Pj{\Ax\^ + \Ay\^). (l7) 

Only the first five quartic terms (/3i — P^) are required to describe '^He The 
additional two terms here (/Jg and Pr) appear because the rotational symmetry 
of the crystal is discrete, where as rotational symmetry is continuous in the 
fluid. Gradient terms can also be calculated [ p^|Jl9| , ^ , but we will not make use 
of these here. 

To second order in A the free energy in a finite magnetic field, Fh, is 
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(18) 



where, 



Q^{k,iu)n) = {iiOn + Sk- IJ. + User ■ H) ^ (19) 



and LUn are the Matsubara frequencies. 
Thus to all orders in H 



(iun +ek - H - HBcr* ■ H) 



Hence, 



where. 



X (cT . A. Sin k.+a.Ay sin ky)a, - ^^.^^ ^ _ _ 
x{a* ■ A*^smkx + (T* ■ A*ySm.ky)a2. (20) 



Fh = A^X K + ^xX A*+Ayx A*+Ayx A* (21) 

=xx =xy =yx —yy 



^S'' = -^E/dfcsinfciSinfc,- (22) 



[{icon - £fe + M)' - \H\^] [(i^n + - M)" - l-ff P] 

By x-y symmetry Xxy = Xyx = 0. Some algebra then leads to 



Fh = [ao + a2\H\^)i\A,\^ + \Ay\^) + ia^H ■ (A, x A*^ + Ay x A;) 

-2a2\H-Af. (23) 



Where, 



^2^/- sin-fc^((£fc-/x)^+0 

"0 pl^J [(i^^ _ + _ Iff |2] [(i^^ + _ _ |ff|2] ' l^^^ 

lUJn 

P [{iun -ek + lif - |2] [{iujn + ek- - \H\^] ^ ' 
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and 



Clearly ag reduces to a ( [T^ ) in zero magnetic field, but the ai and a2 terms do 
not have an analogue in the zero field Ginzburg-Landau expansion. It is inter- 
esting to note the similarity of these extra terms to the change in the Hartree- 
Fock-Gorkov quasiparticle spectrum caused by the magnetic field - (^ij). The 
cross product of any complex vector with its complex conjugate is purely imag- 
inary]^ so the square root of minus one before the ai term in the expression for 
the free energy is to be expected. 

As we have expanded in A but not in H the above expression for the free 
energy is valid for small gaps at all field strengths. It is therefore valid close to 
He- But, note that, since we assumed an exchange-only magnetic field we do not 
consider the vortex lattice here. Agterberg and IIeeb||l^,|2^ have discussed the 
vortex lattice using Ginzburg-Landau theory, but did not include the Zeeman 
terms of (p3|). 

In the Ginzburg-Landau formalism the superconducting phase transition oc- 
curs when the quadratic terms go to zero. In a zero field this condition is simply 

a(r - Tc) = 0. (27) 
In a finite spin only magnetic field the equivalent condition is that the matrix 

£ = ayA,^* (28) 

has (at least) one zero eigenvalue, but no negative eigenvalues, where the indices 
i and j run over both orbital and spin degrees of freedom. In this case 



g= H« ' (29) 




where, 

(ao + a2\H\'^ - 2a2-ffi iai^a -iaiH2 

-iaiHs ao + a2\H\'^ - 2a2Hi 

iaiH2 —\aiHi ao + a2\H\'^ — 2a2H^ ^ 

(30) 

The condition for there being a zero eigenvalue of a is 

^ This can easily be confirmed. Consider the cross product of the most general complex 
vector, V = {a + ib, c + id, e + if). It is trivial to show that v xv* = —2i{cf — de,be — 
af, ad — be). 
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(ao + a2\Hf ~ 2a2Hl) (ao + a2\H\^ - 2a2-ff?) (ao + a2\H\^ - 2a2H^) 
- {ao + a2\Hf - 2a2Hl) ajH^ {ao + a2\Hf - 2a2Hi) a\Hl 
-{ao + a2\H\^ -2a2Hi)alHl = 0. (31) 

This expression can be greatly simplified by choosing our coordinate system so 
that H lies parallel to one of the axes. With, for example, H = (0, 0, H) we find 



iaiiJ 

§_ = \ -\aiH ao + a2i?2 | . (32) 

ao — a2H'^ 




Which has at least one zero eigenvalue when 



(ao - a2i?')((ao + a2Hy - ajH^) = 0. 
The eigenvectors of a are 



(33) 
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1 
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(34) 



Where k is real. To second order in A, k is given by 



^^_ao+^ (35) 
aiH 

Much recent work (see introduction) has suggested that Sr2Ru04 is likely 
to be in an state analogous to the A-phase of ^He. If the pairing interaction 
favours the A-phase in zero magnetic field there are three possible solutions in 
a magnetic field. 



A, = -iAy = (0, 0, 1) (36) 
Ax = -iA,j = {l,iK,Q) (37) 
A,^~iAy^{-iK,l,0) (38) 

Equation^ is the A-phase with d(fe) parallel to H. Equations ^ and ^ both 
give the A2 -phase for < |k| < 1 and the A 1 -phase for |k| — 1. Analogy may 
be drawn to the description of cUiptically polarised light in optics One can 
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think of the three A-hke phases as being described by an elhpse of eccentricity 
Vl — K^. The A-phase is the special case of linear polarization when the ellipse 
reduces to a line parallel to d{k). The Ai phase is the special case of circularly 
polarized light a circle which lies in the 1,2-plane. The A2 corresponds to any 
ellipse between these two extremes. In the Ai and A2-phases by taking the 
appropriate superposition of (|3^) and (|3^) the major axis of the ellipse can be 
made to point in any direction in the plane perpendicular to H . 

4 Numerical results 

To progress further we must resort to solving the self consistent Bogoliubov- 
de Gennes equations numerically. To do this we fit the hopping integral and 
site energy to the experimentally determined Fermi surface of the 7-sheet of 
Sr2Ru04 The interaction potential is restricted to include nearest neighbour 
terms only and chosen to give the experimentally observed critical temperature 
{1.5K). 

4.1 d(fe) parallel to H 

We begin by studying the first solution of the Ginzburg-Landau theory (p6|), in 
which d(fc) is parallel to H. In zero field we find that the ground state of the 
model is a triplet state analogous to the A-phase of "^He, specifically the state is 

d = Z\o(sin/c3; -I- isinfcj,)e. (39) 

Here we have defined the vector order parameter to point in the e direction. 
In zero field, all directions in spin space are degenerate if spin-orbit coupling 
is neglected. When an external field is applied the ground state has d{k) per- 
pendicular to the field, as we will show below. However, in Sr2Ru04 the order 
parameter is thought to be aligned with the c-axis§, by spin-orbit coupling. 
Therefore despite the low critical field along the c-axis, Sr2Ru04 presents us 
with the possibility of studying a triplet superconductor with a magnetic field 
parallel to the order parameter. It is therefore interesting to predict what would 
be observed in such experiments. To do this we simply discard any A-phase 
like solutions with d{k) not parallel to H. We then consider the remaining self 
consistent solution of the BdG equation with the lowest free energy. 

A field applied parallel to d{k) does not cause a change in the symmetry of 
the gap. It follows that at zero temperature the gap is independent of magnetic 
field strength (see appendix). At finite temperature, a field applied parallel to 
the order parameter causes a change in the magnitude of the gap (see Fig. 0.) It 
should be noted that the gap is nodeless but has minima ed — and ky = 0. 

We calculate the heat capacity, magnetisation and magnetic susceptibility as 
functions of temperature and field strength. For an isotropic nodeless gap in zero 
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Fig. 1. (a)Tlie Fermi surface and the gap at T /Tc = 0.5 with (b) fisH/kBTc = 0, (c) 
fiBH/kBTc = 0.5, (d) figH/kBTc = 0.9 

field it is well known that these properties behave as 

a,M,x~exp(-A.). (40) 
Kb J 

We find that for an anisotropic, nodeless, p-wave gap the thermodynamics 
have the same form, even in the presence of a magnetic field (see inset Fig. H) We 
therefore define the effective gap, ^eff 'seen' by the thermodynamic functions 
as 



a,M,x^exp(^). (41) 

We find that ^eff is the mean gap at the Fermi surface, |<i(fcir)| in zero field 
and that ^e// is a linear function of magnetic field strength (see Fig. ^j.) That 
is to say that 

Z\e//=Rfc;^-A^s|Jf|. (42) 
4.2 d{k) perpendicular to H 

Recall that the ground state of the model in zero field is 



d = Z\o(sin kx + i sin ky){l, 0, 0). 



(43) 
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l/T 



0.5 1 1.5 

Fig. 2. Aefj (normalised to |(i(fe_F)| at T = = 0) as a function of magnetic field 
parallel to d{k) extrapolated from heat capacity (circles), magnetisation (squares) and 
magnetic susceptibility (diamonds). The line is |d(fcF)| — hbH. Inset - Logarithmic 
plot of heat capacity with inverse temperature at various fields. From the bottom up: 
H=OT, 0.28T, 0.42T, 0.71T, 0.85T, 1.13T, 1.41T, 1.76T, 2.12T, 2.47T and 2.82T. 

(See Fig. |^a.) We will now examine the numerical solutions of the fuU BdG 
equations corresponding to the second solution of the Ginzburg-Landau theory 
( p7| ) and(^). In a magnetic field the ground state is when the vector order 
parameter points perpendicular to the field. There is also a change in the pairing 
state to a phase analogous to the A2-phase of ^He, where 



d = Zio(sinfca, + isinfcj^)(l,iK,0). (44) 

where k is a real function of temperature and field strength (Fig. ||b,c.) Physically 
this corresponds to the majority of the spin 1 Cooper pairs aligning themselves 
antiparallel to the magnetic field. 

In ^He as the field and temperature increase k increases until k = 1. This is 
the Ai phase which is the ground state of '^He near to Tc in finite fields. The 
A 1 -phase has order parameter 



d = Z\o(sinA:j; + isinfcy)(l,i,0) (45) 

and corresponds to single spin pairing with all of the Cooper pairs aligning 
themselves with the magnetic field (Fig. ^.) However, even near Tc and in 
large fields we do not find that the Ai-phase is the ground state of our model. 
If such a transition does occur then it is certainly well above the experimentally 
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Fig. 3. The Fermi surface (circles), spin up gap (crosses) and the spin down down 
(squares). In the (a) A-phase {T ^ H = 0), (b) the A2-phase (T=0, H=1.4T), (c) the 
A2-phase with a larger k {T — l.SK, H — lAT) and (d) the Aiphase - not observed. 



observed upper critical field. This is in agreement with experiment as no Ai- 
phase has been observed to date. 

Due to the nodeless gap in the A2 phase the specific heat has an exponential 
temperature dependence. Hence we can calculate the effective gap for this field 
orientation (Fig. ^.) We find a linear field dependence in low fields but its de- 
pendence is much weaker than for d{k) parallel to H and there is an upturn in 
large fields. There is known to be a qualitative change in heat capacity in this 
field orientation [|l2j. It remains to be seen if these are related. 

5 Conclusions 

We investigated the order parameter of Sr2Ru04 in an exchange-only magnetic 
field. A Ginzburg-Landau symmetry analysis implied three possibilities: either 
a ^He Al or A2 phase with d{k) perpendicular to the magnetic field or a pure 
^He A phase with d{k) parallel to the magnetic field. We explored the exchange 
field dependence of the order parameter and energy gap in a one-band model of 
Sr2Ru04. The numerical solutions showed no Ai phase for physically reasonable 
field strengths and that of the two remaining phases the A2 phase is lower in free 
energy. We did not include the effect of spin-orbit coupling which could change 
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Fig. 4. Aeff (normalised to \d{kF)\{T = H = 0)) as a. function of H perpendicular to 
d(fe) (solid line) extrapolated from heat capacity. For comparison we plot Af.f / for H 
parallel to d{k) (dashed line). 



the ground state for particular orientations of the magnetic field (particularly 
with H parallel to the c-axis of the crystal.) We investigated the behaviour of 
the heat capacity as a function of both field and temperature for both of these 
solutions. We have shown that the variation of the exponential cutoff below Tc as 
a function of H is quantitatively and qualitatively different for these two phases. 
This makes heat capacity an excellent experimental probe of the symmetry state 
in a magnetic field. 

We acknowledge support of BJP from an EPSRC studentship. 



Appendix 



For an A-phase triplet superconductor with H parallel to d{k) and z the BdG 
equations are 
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/fife + IJiBH 
Ck-usH dsik) 
d^(fe) -e-k-f^sH 



d3{k) \ /u^^{k)\ 
ui„{k) 
v^,{k) 







^E„{k) 



Ula{k) 
V^a{k) 

Wik)J 



.(46) 



Hence, the eigenvalues are 



where 



= Eo{k) + (JHBH 



Eo = Vefe + M3(fc)P 
is the spectrum in zero field. The eigenvectors are 



and 



dsik) 



v/(i?o(fe)-efc)2 + |d3(fc)P 

Eojk) - Ek 

^{Eo{k)~tkY + \d^{k)\^ 



(47) 
(48) 

(49) 
(50) 



Substituting these into the self-consistency condition (^ we find that the gap 
equation is 



d^ik) ^-^U (fe);^ tanh( ^ ) 



fccr 



At T = this becomes 



(51) 



(52) 



h<7 



which is independent of H . 
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